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' Abstract. We derive a three-term asymptotic expansion for the ex- 

I pected lifetime of Brownian motion and for the torsional rigidity on 

^ . thin domains in E,", and a two-term expansion for the maximum (and 

' corresponding maximizer) of the expected lifetime. The approach is 

__ , similar to that which we used previously to study the eigenvalues of the 

■ Dirichlet Laplacian and consists of scaling the domain in one direction 
' and deriving the corresponding asymptotic expansions as the scaling 

. , parameter goes to zero. As in the case of eigenvalues, these expansions 

• may also be used to approximate the exit time for domains where the 

' scaling parameter is not necessarilly close to zero. 

S: 

1. Introduction 

^ I Let be a bounded open set in M" and consider the elliptic equation 

^! -An = 2, xen 

^ ■ (1-1) 

in ■ u = o, x£dn 
in ■ 

I Under certain mild conditions on the regularity of the boundary, the above 

■ equation has one and only one non-negative solution in (7(0). However, and 
0^ . except for a few domains such as ellipsoids, this is not known in closed form 

' - see |AFRj for a recent result in the case of equilateral triangles. 

^> . On the other hand, solutions to equation (jl.ip have a probabilistic in- 

^ I terpretation in terms of the Brownian motion associated to the Laplacian 

^ • in M". More precisely, the value of u at each point x yields the expected 



lifetime of a particle starting from x and it is thus of interest to be able to 
obtain approximations to these solutions. This is the purpose of the present 
paper where we apply to this problem an approach that we have used re- 
cently in the case of Dirichlet eigenvalue problems and which provides quite 
accurate approximations for the first eigenvalue of the Laplace operator for 
certain classes of domains |BFH IBF2j . The idea consists in scaling the do- 
main $7 along one direction and determining the asymptotic expansion of 
the solution in terms of the scaling parameter as it goes to zero. In this 
way we obtain an asymptotic expansion for the solution of equation (jl.ip 
from which it is possible to derive expansions for other quantities such as 
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the maximum of u, the corresponding maximizer and also the integral of u. 
In line with the above interpretation, the second of these quantities corre- 
sponds to the point in the domain with the largest expected lifetime, while 
the last one is known in the elasticity literature as the torsional rigidity - 
see jBBC] and the references therein. As may be seen from the examples in 
Section [H although our approximations are quite accurate for a farily large 
range of values of the scaling parameter, the error can vary a lot as this 
parameter approaches one, depending on the initial domain considered - see 
the discussion in this section. 

It should also be mentioned that there is a large number of papers de- 
voted to the asymptotic expansions of the solutions to elliptic boundary 
value problems in thin domains - see [N], |NTj and the references therein. 
However, and to the best of our knowledge, the problem considered here has 
not been studied from this point of view. 

The paper is organized as follows. In the next section we lay down the 
notation and state the main results of the paper. In Section [3] we prove the 
asymptotic expansion for solutions of the elliptic problem depending on a 
scaling parameter. Sections U] and [5] then present the asymptotics for the 
maximum of this solution and for the torsional rigidity, respectively. We 
finish with an analysis of the error of these approximations for some specific 
domains. 

2. Formulation of the problem and main results 

Let X = {x',Xn), x' = {xi, . . . ,Xn-i) be Cartesian coordinates in R" and 
R"~^, respectively, a; be a bounded domain in 1R'"~^ with C^-boundary. By 
h± = h±{x') S C{uj) n C^{oj) we denote two arbitrary functions such that 
H{x') := h-{x') + hj^{x') > for x' G a; and H{x') = on du. We also 
define the two functions 

d{x') = h+{x') - h^{x') and p{x') = h+{x')h_{x') 

We introduce a thin domain by 

Vie := {x : x' £ to, — e/i_(x') < Xn < £h^{x')}, 

where e is a small parameter, and consider the problem 

(2.1) - Aue = 2 in Q^, Ue = on dQe- 

In view of the smoothness of the functions h± and the boundary of uj the 
domain satisfies the exterior sphere condition at every boundary point. 
By Theorem 6.13 in \GT\ Ch. 6, Sec. 6.3] it implies that the solution to the 
problem (|2.ip belongs to (^(Og). 

Let '■= Xn/s- Assuming that the functions h± are smooth enough, we 
introduce a sequence of functions 

[/2(x',e):=4'^(x')d, 
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^2 3^ «r(^'):=P(^'), a['\x'):=d{x'), 



af\x') := -1, 



(2.4) U2j{x',^n) := U2j{x',U) + + a'^'Hx'), 

2i-i 

(2.5) U2j{x',C,):=Y,c^?'\^')C, 

i=2 

(2.6) af^)(x') := --^^A^^a'fS'Hx'), i > 2, 

(2.7) af^\x') := - ^ ap)(xO {Ki-'Ti " h.{x')Y—\ 

1=2 m=0 

(2.8) a^o^\x') := -U2j[x' ,h+{x')) - af^\x')h+{x'). 



The first of our results describes the uniform asymptotic expansion for 
and forms tlie basis for the remaining expansions in the paper. 



Theorem 1. Given any N ^ 1, let the functions h± he such that 
Q/p-') G C^(lIJ), j ^ N . Then the function satisfies the asymptotic formula 



(2.9) u,{x) = £e^^'^2, (x', ^) + 0(^2^+2) 



in the C{^}i;)-norm. In particular, 
(2.10) 

M^',^n) = -^l + ^nd{x')+p{x'), 

Mx',(n) = -§A,,d{x')-§A,,p{x') + a^i\x')^n + al^\ 
Mx',(n) = ^Al,d{x') + ^Al,p{x') 



||A,, { [d^x') +pix')] A,,dix') + 3dix')AMx')} 

-f^A^. [d{x')pix')A^>d{x') + 3p(x')A^/p(x')] 
+af)(x')^„ + 4')(x'), 
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where 

a^^\x') = l[£{x')+p{x')]A,,d{x') + ^d{x')AM^'), 

al^\x') = ldix')p{x')A,,d{x') + ^pix')AM^'), 

af\x') = ^d{x')A^,[d{x')p{x')A^>d{x')+3p{x')AMx')] 

+^ [dHx')+p{x')] A,, {'id{x')AMx') + [d^x') +p{x')\ A^.d{x')] 
-^d(x') + 2p{x')\ AIp{x') 

[d\x') + M^{x')p{x') +p\x')\ Al,d{x'), 

o^t\^') = ^^.'[d{x')p{x')A^,d{x') + ?,p{x')AMx')] 

^ d(xW) ^^,|[^2(^/) +^(^/)] A.,,d{x')+M{x')AM^')] 
-^p{x') [d\x')+p{x')\ Al,p{x') 
-^d{x')p{x') [d^{x') + 2p{x')\ Al,d{x'). 



The coefficients of the asymptotic expansion (j2.9p involve two scales, 
namely, the variable x' and the rescaled variable so, this is a two-scale 
asymptotics. This is a very natm'al situation for the problem (j2.ip since by 
passing to the variable we get a bounded domain 

^ ■■= {{x',(n) : x' E w, -h-{x') <^n< h+{x')}. 

In this domain there isn't a distinguished variable as was the case of x„ in 
the domain O^. This is one reason why the asymptotics for involve two 
scales. Another way of understanding this fact is that while Xn ranges in 
a small interval the remaining variable x' ranges in a finite set. From this 
point of view, it is natural to rescale the variable Xn and pass to 

From Theorem [H we are able to derive explicit asymptotic formulas for 
both the maximum of and the torsional rigidity for the family of do- 
mains as e goes to zero. 

Theorem 2. Let the family of domains fig and the functions /i+ and H 
be as above. We assume further that H satisfies the following hypothesis 

HI There exists a unique point x £ uj at which H achieves its global 

maximum which will be denoted by Hq; 
H2 The Hessian matrix of H at x, denoted by 2H2, is negative definite; 
H3 The functions h± are 5 times continuously differentiable in a vicinity 

ofx. 

Then the maximum value of in for sufficiently small e satisfies 

max Ue{x) = + \hI [4 tr + 2 tr P2] + 0{e^), 

xefij 4 o 
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as e 0. Here do is the value of d at x, while the matrices 2D2 and 2P2 
are the Hessian matrices of d and p at the point x, respectively. 

Remark 2.1. If one drops the hypothesis that H2 is nonsingular it wih still 
be possible to obtain an expansion, but this will be much more involved and 
will depend on higher order terms in the expansion of H around x. 

Remark 2.2. The hypothesis on a unique maximum of H may also be 
dropped and provided there is only a finite number of such maxima the 
results still hold except that one has to construct different expansions for 
each maximum. The case of a continuum of maxima is also possible to han- 
dle with the techniques employed here but requires some further changes to 
the approach. 

Remark 2.3. In the process of proving the above theorem we also obtain a 
two-term asymptotic expansion for the maximizer - see Theorem |H How- 
ever, this depends on higher order terms in the expansions of both d and p 
and requires the introduction of more detailed notation which we postpone 
till Section [4] below. 

Remark 2.4. Under the hypotheses HI and H2 of Theorem [2] and assuming 
that the functions h± are smooth enough in a vicinity of x, it is possible to 
construct more terms in the asymptotic expansions for the maximum of in 
and for the corresponding maximizer. In order to do it, one should follow 
the main lines of the proof of Theorem [2l employing Lemma [4. 2 1 as a starting 
point. At the same time, it requires bulky and technical calculations which 
we would like to avoid. This is the reason why we provide only two-terms 
asymptotics in Theorem [2j 

Finally, the integration of the asymptotic expansion for given by Theo- 
rem [1] yields the corresponding asymptotic expansion for the torsional rigid- 



Theorem 3 (Torsional rigidity). Under the conditions of Theorem[l\ we 

have 



ity. 




+ {HHy)d\y) +p{y)H{y) [p{y) - d\y)\ } A2p(y) 

- ^^[H\y)d{y) -2d{y)p{y)H{y)\ 
X Ay { [d^{y) +p{y)] Ayd{y) + 3d{y)Ayp{y)} 
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-j^[HHy)-My)H{y)] 

X Ay [d{y)p{y)Ayd{y) + 3p{y)Ayp{y)] 
+ y{y)H{y)af\y) + Hiy)a^^\y) dy 

where a^^^ and a^^^ are as in Theorem\^ 

Remark 2.5. In order to obtain the average expected lifetime it remains to 
divide by the volume of Q.^ which is given by 



/ dx= dindx' = e H{y) dy. 

JCl^ JujJ—eh—(x') J u) 



3. The asymptotic expansion for 

In this section we prove Theorem [TJ We begin by passing to the variables 
(x', in (j2.ip leading us to 

(3.1) (^-e^A^/ - Ue = 2e^ in = on dVL. 
We construct the asymptotic expansion to the problem ()3.ip as follows 

oo 

(3.2) Ue{x) = Y,e''^U2j{x',in), 

j=0 

where Uj{x',(,n) are functions to be determined. 

We substitute the expansion (|3.2p into (j3.ip and equate the coefficients of 
like powers in e. This yields the following boundary value problems for U2j: 

,3.3) -W^'- 

^^2 = 0, ^n = ±h±{x'), x'GoJ, 



d' 



U2j 



U2j=0, Cn = ±h±{x'), x'guj, j ^ 2. 



It is easy to check that the solution to the problem (j3.3p is given by the 
formulas (12. 2p . (12. 3p . Given U2j-2, the solution to the problem (j3.4p is 
determined uniquely. This is why we can find the functions U2j in a recurrent 
way. Moreover, at each step the solution to the problem (13. 4p can be found 
explicitly. These explicit solutions are given by formulas (12. 4p . (j2.5p . (j2.6|) . 
(|2.7p and (|2.8p that can be checked by direct calculations. 
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Given any N 1, assume that a^^"'^ G C^(aJ), j ^ N. Let 

N 

(3.5) u^{x',^n):=Y.e'^U2j{x',^n)- 

3=0 

It follows from the problems (j3.3p . (j3.4p that the function solves the 
boundary value problem 

(^-e'^^x' - = + e^^+^A^/naTv in n, = on d^. 

Hence, the function := — is the solution to 

(3.6) (^e^A^> + = e^^+^A^,U2N in = on dQ. 

The coefficient affecting the derivative in the last equation is one. Em- 
ploying this fact and applying the maximum principle in the form of in- 
equality (1.9) in [LIU Ch. 3, Sec. 1], we obtain the estimate 

(3.7) ||5f Ib(iT) ^ Ce''''-'\\A.'U2N\\c(n) ^ Ce^''^^ 



where the constant C is independent of e. It proves the formula (|2.9p . The 
formulas (f2T0]) follow directly from ([23]), (f277l) . 

(j2.8p . The proof is complete. 

4. Proof of Theorem [2] 

In the whole of this section we shall consider the function Uf,{x' , ^n)- This 
is then defined on and it is clear that it is sufficient to find the maximum 
of «e(x',i^„) since after rescalling x„ — > x^e the maximum of the function 
remains unaltered. 

4.1. Existence of an expansion and terms of order e^. We begin by 
showing that, under the hypothesis of Theorem [2] and up to order e^, the 
maximum of u has an asymptotic expansion that may be obtained directly 
from the expression of U2- 

Lemma 4.1. Under the hypothesis of Theorem\M we have 
maxne(x) = maxu2{x)e^ + ©(e^) as e — > 0. 

Furthermore 

{_ d{x) \ 1 2 
maxn2(x)=.2 ^ J = 4^0 

and is unique. 

Proof. The first part follows directly from the asymptotics of given by 
Theorem [H since we now have 

Ue{x',Cn)=e\2{x',U) + 0{e^). 
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For the second part, note that we may write 

U2{x\in) = -[en + in{h_{x')-h+{x')) + h_{x')h+{x')\ 

= - in-\d{x)\\\H^{x). 

This last expression is clearly maximized when ^„ = d{x)/2 and H is also 
maximized, yielding x' = x and ^ = d{x)/2. The uniqueness follows directly 
from hypothesis HI. □ 

In order to go on to obtain the next terms in the expansion for the maxi- 
mum (and the corresponding maximizer), we need to show the existence of 
such an expansion which we do in the next lemma. This also proves that 
the coefficients of the terms of order one in both the expansions for x' and 
^„ vanish. 

Lemma 4.2. Given any N ^ 1, assume that the functions h± are [A^/2]-|-2A^ 
times continuously differentiable in a vicinity of the point x, and the hypothe- 
ses HI and H2 of Theorem\^hold true. Then the function u^^n has only one 
stationary point which is a maximum. The corresponding maximizer has the 
following asymptotic expansion 

[N/2] 

]\T — X 



(4.1) 



e,N 

i=l 



[N/2] 
i=l 

The maximum of the function Ug satisfies the identity 

(4.2) maxne(a;) = maxue(a;', e^„) = maxnf (x', ^„) + 0(5^^+^), 

and any maximizer (x^,^^) of this function has the asymptotic expansion 

[N/2] 



(4.3) 

W/2] 



i=l 



Let us find the maximum of . In order to do this, we should first find 



Proof. The identity (|4.2|) follows directly from the asymptotics (|2.9|) for u, 

n of . In order to d 
is functions by solving 

V(.',5„)uf(x',en)=0, 



the stationary points of this functions by solving the equation 



which is equivalent to 

N 

(4.4) 5^e2^-2^(,,^„)tX2^(x',^„) = 0. 



i=l 
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It follows from Lemma HTT] that for e = this equation has a unique solution 
(x, d(x)/2). In order to solve (j4.4p for e > we apply the implicit function 
theorem considering as functions of e. We first need to check that 

the corresponding Jacobian is non-zero. It is easy to see that this Jacobian 
at £ equal to zero coincides with the determinant of 2H2 which is non-zero 
by hypothesis 112. 

The assumption for the smoothness of h± and the formulas (12. 2p . (j2.3p . 
(1231), dlS]), (I2S1), (1221), for U2i yield that U2i{x',in), i = l,...,N, are 
[A^/2] + 2 times continuously differentiable in a small vicinity of (x, d{x)/2). 
The dependence of the left hand side of (j4.4p on is holomorphic and by the 
implicit function theorem we conclude that for e small enough there exists 
a unique solution (x^ ^, £,s,n) to (|4.4p which is [A^/2] + 1 times continuously 
differentiable in e^. Hence, we have the Taylor polynomial (j4.ip . The point 
(x^ j'^t^cn) is the maximizer for since by hypothesis H2 the Hessian of 
at this point differs from 2H2 by an error of order 0{e'^). We employ 
this fact and expand (x',^„) in Taylor series at (x^^,^e^7v)- As a result, 
we have the estimate 

(4.5) nf (x', en) - (4,7V, Ce,N) ^ "Cl (|x' - x',^^\^ + j^n " ^..^vl') , 

where Ci is a positive constant independent of e, x' and ^n- This estimate 
is valid in a small fixed neighborhood Q of the point {x,d(x)/2). Since the 
function has the maximum at (x^ N,S,e,N), we can choose the neighbor- 
hood Q so that outside it the estimate 

(4.6) (x', Cn) - «,7V, ^e,N) ^ -C2 < 

holds true, where the constant C2 is independent of e. Let us choose (x', (,n) 
so that 

(4.7) |^'_4_^|2 + |^^_^^_^|2^C73e^A.+2^ 

where C3 is a positive constant independent of e, x' and Then it follows 
from (14. 5p . (14. 6p that for such (x^^n) the inequality 

nf (x',e„) nf (x^,^,ee,7v) - CiC3e2^+^ 

is valid. Together with (j4.2p it implies that a maximizer (x^, ^e) of can not 
satisfy (j4.7p for sufficiently small e and sufficiently large C3 and therefore 

(4.8) |4_4_^|2 + |^^_^^_^|2^C352A^+2_ 

This inequality and ()4.ip prove (j4.3p . □ 

4.2. The terms of order e^. In order to determine X2 and ^2 we shall need 
the terms of order in the asymptotics of the gradient of , for which we 
need to consider U4. We must also develop d and p around x. In full 
generality, and to obtain the full asymptotic expansion, these developments 
should be written in terms of homogeneous polynomials of increasing degree. 
However, and since in order to obtain the first two terms in the asymptotics 
we will only need terms up to the homogeneous polynomials of third degree. 
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we shall choose a form that will be more convenient for our calculations. 

Write thus d and p as follows. 

(4.9) 

d{x') =do + d{{x' -x) + (x' - xYD2{x' -x) + D^{x' -x) + 0{\x' - x]^), 
p{x') = po+p\{x' -x) + {x' -xfP2{x' -x) + Pzix' -x) + Odx' - x|^), 

where di = Vx'd{x), pi = V^'pix). 

Due to the relation between d and p via the functions h± and the fact 
that H'(x) = h'_{x) + h'j^{x) must vanish, we easily obtain that 

Pi = - ^dodi. 

The homogeneous polynomials of degree three, and P3, will be specified 
below. 

In the case of U4, the relevant derivatives are given by 

V.U4(x',^„) = - le^V,, [A,,(i(x')] - ^e^v^.' [A,,p(xO] 

HnV.>a^^\x')+Vx,a^^\x'). 

We shall first obtain the term of order in the derivative of with respect 
to This will have a component coming from the term of order in the 
corresponding derivative of U2, and another from the constant term in the 
derivative of U4. In the first case it is straightforward to obtain that the 
required coefficient is given by 

(4.10) - 2U + d\x2. 

In the case of U4 the term coming from a[^^ is given by 

a^^\x + 0{e^)) = ^ [dl+po) ii{D2) + dotr(P2) + 0{e^). 
We thus obtain 

|g(x + 0(e2), 14 + 0(^2)) = -\dltT{D2) - dQiT{P2) 



+\ {dl+po) tv{D2) 
+4tr(P2) + 0(e2) 
^ {dl + 4po) ir{D2) + 0{e^) 



This, together with (j4.10p . yields 

|^(x + X2e^ + 0(e3), 1 do + + 0(e=^)) 
(4.11) "^^^ 



-2U + d\x2 + ^{dl + 4po) tv{D2) + 0{e^). 
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We will now proceed to compute the gradient with respect to x' . The 
case of U2 is again straightforward and we obtain 

(4.12) = (^1^0 + ^„2e2) [di + 2D2X'ie^\ - \d^dx + 2P2X2e'^ + 0{c^) 
= {iradx + d^B2X2 + 2P2X2) + 0(^3). 

In the case of we are only interested in the terms of order e*^. However, 
there are now expressions of the form 

Vx.'A^/d(a;') and V^.'A^./p(x'), 

this being the reason why we need the homogeneous polynomials of third 
degree in the expansions of d and p. On the other hand, this implies that 
the only relevant terms from D3 and P3 are those where one of the variables 
appears at least twice. If we write 



D^{x') 



n-1 

^ ^ dijkXiXjXf^ 
ijk 

n-1 



E 

i=l 



n-1 



diiiX^ 



with 



+ ^ ^ {diij -\- diji -\- djiijXj^Xj 

n-1 



+ ri{x'\ 



''3 ("^ ) — ^ ^ dijf^XiXjXjf, 

we may assume without loss of generality that the coefficients dijk are in- 
variant under any possible permutation of the indices. If we then denote dm 
and diij = diji = dja {i 7^ j) by Sa and respectively, the expression for 
Ds becomes 



n-1 



n-1 



j=i 



With this notation we get 

V.,,{^^,D^{x')) 



6V, 



n-1 



n-1 



i=l 



'n—1 n—1 



12 
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In a similar fashion, if we write 



n-1 



^^{x') = ^ j TTiixf + 3 ^ -^ijXiXj j +rf(x'), 



1=1 



we get 



n-l 



n-1 



i=l 



'n-1 



n-1 



In this way, we obtain after some lengthy but straightforward calculations. 



do 
T 



do 



+ 



do 



(dg +Po)S5 + ^^{2dodi - ^dodi) + 3doS^ + trPac^i 



+doPoS5 + ^^{-\dldi +podi)+3po^n - ^dodi + 0(e2) 



trPo 



If ^2 
3 



{dl + 4po) 



+ + 



where 

^5- 



'n-1 



n-1 



'n-1 



trDo 



n-1 



^ (^ji, ■■■,Y1 '^J>-i and = ^ ^ji, . . . , ^ TTj,n-i ■ 



Combining this with (j4.12p yields 



(4.13) 



in2dl + doD2X2 + 2P2X2 

+ l{dl + Apo) [3S^ + ^0^5 + ^c?i 



from which we obtain the second equation for X2 and by equating the 
coefficient of to zero. From equation (j4.1ip we get 

(4.14) en2 = ld\x2 + ^{dl + 4po) tr D2. 

Substituting this into equation (j4.13p yields 

^ {didi) + 2P2 + doD2]x2 = -^(d2 + 4po) (^a!i + 3S, + §doS5). 

To prove that there is a unique solution, we must show that the matrix 
multiplying X2 is nonsingular. In order to do this, we shall relate the terms 
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apearing above to those in the expansions of the functions h± and H. If we 
write 

h±{x') = h^ + (/if)* (x' -x) + {x' - xfHf{x' -x) + 0{\x' 
we see that 

do = /iq ~ K 
di = ht — /i7 

D2 

H2 



H. 



2 -'-'2 

+ Hr, 



and 



Po 



pi = h^h^ + hQ hi 

P2 = h^H^ + h^H^ + ht{h^y. 
Replacing this in the expression above yields, after some manipulation 

doD2 + 2P2 + ^did\ = HoH2 + ^{ht -h^){ht -h^Y + htiKY 

= HqH2, 

where we used the fact that = Hi = hf + . Since the matrix HqH2 is 
negative definite by hypothesis, we may invert it to obtain 



(4.15) 



X2 



1 



trD 



-!2-^H2 ^di + 3i?2 ^ ('^TT + ■^d-o'^s 



where we have used the fact that + Apo = Hq. Replacing this back 
into (|4.14p yields 



(4.16) 



Cn2 



^ tiD2d\H-^di + 34^2"^ (^n + ^do^. 



+^HitrD2. 



If we now evaluate U2 and U4 at the maximizer we obtain 



and 



U2 [x + X2e^ + O{e^),^do + + 0{e^] 



U4 (x + X2e'^ + 0{£^),-do + U2e'^ + 0{e^] 



= -Hi [dotr(D2) + 2tr(P2)] + 0{e^), 

where ^n2 and X2 are given as above. We have thus proven the following 

Theorem 4. Under the conditions of Theorem\^ the maximizer {x'^^,^n*) 
of given by Lemma satisfies the asymptotic expansion 



{x'^^in*) = (^ + X2e^ ^do + Cn2e^) + 0{£^) as £ 



0, 
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where X2 and ^^2 are given by j^. and ( |^.i6p , respectively. The corre- 
sponding maximum satisfies the asymptotic expansion given in Theorem\^ 



5. Proof of Theorem [3] 

To prove Theorem [3] we need only to integrate in fig. More precisely, 
and using the expression given by Theorem [H we have to compute 



/ Ue{x) dx = 
Jq^ J lu J —eh^(x' 



eh^(x') V e / V e 



£3 



+£^uq [x\ ^) + 0(6^) dx„dx' 

I I U2{x\in) dindx' 

Jw J-h^(x') 



+6^ 



i^(x') 

rh+{x') 



r rn+(x ) 

/ / Ui{x',(,n) d£,ndx' 

Jul J~h-{x') 

Jw J~h^(x') 



We shall now compute these three integrals separately. Using the expression 
for U2 we have 



i-h+ix') rh+(x') 

/ U2{x',U)d^n = / -Q + ^nd{x')+p{x') dU 

J-h^(x') J~h^{x') 



11 1 h+ (x') 

lc3^ ^c2^(x') + v(x')£„ 



= -^en + ^end{x')+p{x')^„ 

= -^[hl{x') + h^^{x')] 

+ ^ [hl{x') - hl{x')] d{x') 
+ [h+{x') + h^{x')]p{x') 

and, taking into account the expressions for both d and p, we get that the 
above equals 

I [hl{x') + h'Lix') + 3hlix')h.ix') + 3h+ix')hlix')] = \h\x') 
b D 



as desired. 
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For U4 we have 



h-{x') 



h+(x') 1 1 

-en^x'd{x') - -C^x'Pix') 



h-{x') 



6 



+Q(')(x')Cn + 4')(x') de„ 
-^e„AMx')-len^x'p{x') 



(4)/ 



h+{x') 
-h-{x') 



Using the expressions for Og^^ and a^^ we see that the above integral equals 
_^ [hXix') - ht[x')\ A,>d{x') - 1 [hlix') + hl{x')] AM^') 
+ ^ [{d^x') + p{x'))A^,d{x') + 3d{x')A^,p{x')] [hlix') - hl{x')] 
+ 1 [d{x')pix')A,,d{x') + 3p{x')AM^')] IKix') + /i-(x')] 
= ^ [hXix') - ht{x') + 2hl{x')h^{x') - 2h+{x')hl{x')\ A^,d{x') 

+ ^ [h%{x') + /ii(x') + 3h\{x')h-{x') + 3h+{x')hl{x')] A^,p{x') 
= ^ [d3(x')i/(x') + 4hl{x')h-{x') - 4h+{x')hl{x')] A^,d{x') 

^H'{x')AM^') 
= ^d{x')H{x') [d'^ix') + 4p{x')] A^,d{x') + ^H'\x')A^,p{x') 
= ^H^{x') [d{x')A,,d{x') + 2AMx')] . 

Although the expression for uq is much more involved, the necessary com- 
putations are similar to those above. We need to compute 



h+{x') 
-h-(x') 



UQ{x',^n) d^n 



h+{x') t4 

^Al,d{x') + ^Al,p{x') 



h^[x') 



120 



24 



§|A,. { [d^x') + pix')] A,,d{x') + 2,d{x')A^,p{x')} 



^2 

-f§A^, [d{x')p{x')A.^,d{x') + 3p{x')A^,p{x')] 

+af\x')ir, + Ji\x') de„. 
Using now the following identities 

h%{x') + h^_{x') = H^{x') - 3p{x')H{x') 

hlix') - ht{x') = H^{x')d{x') - 2d{x')p{x')H{x') 

h%{x') + ht{x') = H^{x')d^{x') + p{x')H{x') [p{x') - d^{x')] 

h%{x') - ht{x') = H^{x')d^{x') + 3d(x'y (x')-H'(x'), 
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after some computations we obtain 

rh+ix') 

h_(x') '^^ 

+3d{x')p'^{x')H{x')] Al,d{x') 
+^ {H'{x')d\x') 

+p{x')H{x') [p{x') - d^{x')] ] Al,p{x') 

-TM [HHx')d{x') - 2d{x')p{x')H{x')] 

xA^>{[d^{x') +p{x')] A^>d{x') + 3d{x')A^>p{x')} 

-^[H%x')-3p{x')H{x')] 

xA^./ [d{x')p{x')A^,d{x') + 3p{x')A^,p{x')] 

+ldix')H{x')af\x') + H{x')ai^\x'). 

This, together with the expressions for the integrals of U2 and U4 given 
above yields the formula in Theorem [3l 

6. Examples 

Let us now apply our results to some special cases in order to test the 
accuracy of the approximations for concrete examples. As may be seen from 
Figures 16.11 and 16.21 below, although the error for either the — norm or 
for the torsion stays below 5% for e up to 0.6, this can vary substantially as 
the scaling parameter approaches one. For the examples considered below 
the error at e equal to one for the torsion, for instance, varies between less 
than 2% and 100% in the cases of the folium and the disc, respectively. The 
reason for this is simply that, even in the case where a Taylor (or Laurent) 
series exists for the quantities under consideration, the series expansion for 
these quantities will have a specific radius of convergence. In the case of the 
disc, for instance, we have that the torsion is given by 

vr 



2 1^-2 



e 

which has a radius of convergence of one, thus explaining the large error 
found in this case. 

6.1. Descartes's folium. We consider the case of the domain defined by 
= {{x,y) G m2 . x(x2 + 3e"2y2) _ 3.2 _^ £-2^2 ^ < x < 1} , 

where we have chosen coordinates it in such a way that the scaling is done 
along the y— axis and the corresponding height with respect to this axis is 
minimal. In this case we then have 

H{x)=2h+ix) = 2h_ix)=2xJJ^, dix) = 0, p{x) = ^\~'^f , 
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yielding 

e 
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(1 — x)x^ 
l + 3x 



(6.1) + 



[(3x + 1)=^ - 4] [jx - + (3x + 1)^^ 



[( 



X 



3(1 + 3x*) 
l)x2 + (3x + l)^2^Pi(x^^6 + o(e8) 



(3x + 1)'^ 



where 

piix, := 1 - 30x + 51x2 _^ 4g^3 _^ ^35^4 _^ ^g2x^ + 81x^ - 12^^ _ ^q^^2_ 

The direct apphcation of Theorem El where we only have the explicit 
formula for terms up to order e^, yields 

maxue(x,y) = -(2V3-3)e'^ + -(12 - 7V3)e^ + O(e^). 
xsHe 9 9 



Note that in this case determining the maximum directly from (|6.1|) in order 
to obtain a better approximation implies solving an algebraic equation of 
degree nine. Actually, even solving explicitly for the maximizer using u up 
to order and taking into consideration that we know beforehand using 
symmetry that must vanish, implies having to solve an algebraic equation 
of degree five. 

Computing the expansion for the torsion will in turn yield 



(6.2) 



u{x, ^)(ix(i^ 



f 167r _ 1 U3 
V243^/3 ^' 



\243V3 



+ 



f SOtt 593 \ ^7 



In Figure [6TT] we show, for various values of e, the relative errors for the 
norm and for the torsion of the diference between the values of numerical 
solution determined using the method of fundamental solutions (MFS) and 
the asymptotic expansions given by (|6.ip and (|6.2p . We note that the error 
at £ equal to one - which corresponds to the actual folium - is of the order 
of 3.5% and 2% respectively. 

6.2. Lemniscate. We consider the domain whose boundary is the lem- 
niscate defined by 



X +e y ] 



~2 2 

£ y , 



x~ — £ ~y~, X > 0. 
The functions H, h±, d and p are given by the formulas 

"1 o 1 



H{x) = 2h+{x) = 2/i_(x 
d{x) = 0, p -- 



X2 + 



(1 + 8x2)1/2 



1/2 



i-x2 + i(l + 8xY^ 



UJ 



(0,1) 
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Figure 1. Graphs of the relative errors for the norm 
and the torsion in the case of the folium; the comparison is 
between the numerical solution obtained by the MFS and the 
asymptotic expansion given by equations (16. ip and (16. 2p . 



Some straightforward calculations give 

[(1 + 8x2)3/2 - 2] (l + 2x2 - Vl + 8x2 + 2^2 



+ 



2(1 + 8x2)3/2 



32x^ 



(1 + 8x2)7/2 



+ e 



2 (-512x'5 - 192x^ - 216x2 + 7)\/l + 8x2 _^ 256x'' + 328x2 - 8 



^1 + 2x2 - ^l + 8x2)p2(x) ' 



2(1 + 8x2)7/2 



(1 + 8x2)7/2 

|e'^ + 0(e«), 



P2(x) := \/l + 8x2(512x^ + 192x^ + 152x2 - 5) - 128x^ - 268x2 + 6. 
The maximum of H is now situated at x = ^/2>/ {2^/2) and using 

3 

Ho = ^, do = trD2 = 0, 2trP2 = -2, 
in Theorem [2] then yields 

^2 3e4 
maxue(x) = — - — - + 0{e^). 
xeUe 8 32 

By applying Theorem[3]we arrive at the asymptotics for the torsional rigidity 



Ue{x) dx 



37r - 8 



48 



v3 , , /—^ 37r 
-ln(2 + v/3)-- 



+ 



13 57r 20^3, 

— + — In 2 + V3 

18 16 27 ^ ^ 



e' + 0(e'), 
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where the integrals appearing in the coefficients can be calculated by the 
Euler substitution 



2t Vl + 8x2-1 



vT+sx^ = + 1, x = ^, t 

8 — X 
As in the previous example, we show in Figure 16.21 the relative errors 
for the norm and the torsion in this case. However, comparing the two 
examples gives that for e larger than approximately 0.4 the errors become 
much larger than in the case of the folium. 



Figure 2 . Graphs of the relative errors for the norm and 
the torsion in the case of the lemniscate, with the different 
quantities computed in a similar way to what was done in 
the folium case. 



6.3. Ellipsoids. As mentioned in the Introduction, these are one of the few 
examples where the explicit solution of the corresponding equation (jl.ip is 
known. More precisely, if we consider ellipsoids defined by 



E 



X G 



+ 



+ 



ar, 



1 



we have that the solution of equation (jl.ip in this case is given by 

2 / \ 2' 

i I Xl \ 

U{X) 



1 



1 



Ol an 



1 



an 



The maximum is thus localized at the origin and is given by 

1 



M 



1 



+ ••• + 



1 



Let's assume that a„ is the smallest of the a^'s, and we thus pick this direc- 
tion to be that along which we scale the domain. We then have 

^.,1/2 

iJ(x') = 2/i+(x') = 2/i_(x') = 2o„ 1 



Xn—l 



an-1 , 
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while d vanishes and 



p{x') 



1- — 



^n—l 



From this it fohows that the maximizer is 0{e^), while the maximum has 
the expansion 

n-l 

M 



1 



A straightforward analysis of the error 



E 



M 



n-l 



1 



2 2 I 4 

^ I 



+ 



1 «l 



- + 



+ 



yields that this satisfies 



< < e^(n-\) 



with equality on the right-hand side being achieved for the ball. 

For the sake of comparison with the previous two-dimensional examples, 
we shall now consider the case of the planar disc (in the above notation, n = 
2,ai = 1,02 = e)- In this case the above expressions yield that the relative 
error of the I? norm and of the torsion are, respectively, e^^ -|- 0(6^^) and 

-|- 0[£^\ Thus, although the approximation is very good for sufficiently 
small values of e, the error does become quite large as e reaches one. 
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